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Abstract 

Andreev scattering of thermal excitations is a powerful tool for studying quantized vortices 
and turbulence in superfluid ^He-B at very low temperatures. We write Hamilton's equations 
for a quasiparticle in the presence of a vortex line, determine its trajectory, and find under wich 
conditions it is Andreev reflected. To make contact with experiments, we generalize our results to 
the Onsager vortex gas, and find values of the intervortex spacing in agreement with less rigorous 
estimates. 

PACS numbers: 

67.40.Vs Quantum fluids: vortices and turbulence, 
67. 30. em Excitations in He3 
67.30.hb Hydrodynamics in He3 
67. 30. he Vortices in He3 



I. MOTIVATION 



Superfluid turbulence consists of a disordered tangle of quantized vortex filaments which 
move under the velocity field of each other [l, 2]. If the temperature, T is sufficiently smaller 
than the critical temperature, T^, then the normal fluid can be neglected and the vortices 
do not experience any friction effect [sl. The simplicity of the vortex structures (discrete 
vortex lines) and the absence of dissipation mechanisms, such as friction and viscosity, make 
superfluid turbulence a remarkable fluid system, particularly when compared to turbulence 
in ordinary fluids. Current experimental, theoretical and numerical investigations attempt 
to determine the similarities and the dissimilarities between superfluid turbulence and or- 
dinary turbulence. Questions which are currently addressed concern (i) the existence of a 



Kolmogorov energy cascade at length scales larger than the typical intervortex spacing 



B,3, 



(ii) the existence of a Kelvin wave cascade at length scales smaller than the Kolmogorov 
length 0, 0, y, 3 

(iv) the possible existence of a bottleneck jl2l. Q between the Kolmogorov cascade and 
the Kelvin wave cascade, (v) the nature of the fluctuations of the observed vortex line 
density 1^, [isl, Q, and (vi) their decay 18, [l9], (vii) whether there are two forms of 
turbulence [2^, a structured one, which consists of many length scales (Kolmogorov tur- 
bulence), and an unstructured, more random one (Vinen turbulence), (viii) the effects of 



rotation on turbulence 
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24l |. Most of these questions refer to the important limit 



T/Tc ^ 1, where fundamental distinctions between a perfect Euler fluid and a superfluid 
becomes apparent [25 1. 



Superfluid turbulence experiments are currently performed in both ^He 



and in ^He-B 
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30|]. In the last few years it has been recognized that, to make 



progress in answering the above questions, it is necessary to develop better measurement 
techniques which are suitable for turbulence in quantum fluids. In ^He, the application of 
the classical PIV method 
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331] was a breakthrough. In ^He-B, a non-classical, powerful 



measurement technique which is suitable in the limit T/Tc ^ 1 is the Andreev scattering 29l |. 
developed at the University of Lancaster. 

This article is concerned with the Andreev scattering. The plan of the paper is the 
following. In Section II we shall describe the basic ideas behind the Andreev scattering and 
review what is a quantized vortex line. In Section III we shall write down the governing 
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equations of motion. In Section IV we shall determine the ballistic trajectories of excitations 
in the vicinity of the velocity field of a vortex line, and, in Section V, we shall study the 
transport of heat by ballistic quasiparticles through a tangle of vortices. Section VI will 
apply our result to the current experiments. Finally, in Section VII, we shall draw the 
conclusions. 



II. ANDREEV SCATTERING AND QUANTIZED VORTICES 

The study of the motion of quasiparticle excitations in a superfiuid was pioneered by 



AndreevSJ]. Consider an excitation which moves in the direction of increasing excitation 
gap. The excitation propagates at constant energy, and gradually reaches the minimum of 
the rising excitation spectrum, where its group velocity becomes zero. Thereafter it retraces 
its path but as an excitation on the other side of the minimum. An incoming quasiparticle 
is thus reflected as a quasihole and an incoming quasihole is reflected as a quasiparticle. The 
effect is a consequence of the fact that the minimum of the energy spectrum of the excitation 
lies at nonzero momentum. 

The case of p-wave triplet pairing appropriate to superfiuid ^He has been discussed by 
various authors studying the interaction of excitations with the boundaries [ssl], motion of 
quasiparticles through the A-B phase boundary in ^He[3(|], ballistic motion of quasi-particle 



in slow varying textural field of He-A|37], scattering of ballistic quasiparticles in He-B by 



a moving solid surface [38|, l39|, and calculation of the friction force on quantized vortices [401. 



4l|. Ref. 4l| and [42] are concerned with Andreev reflection within the vortex core and 
therefore apply to the bound states. Our concern is the propagation of thermal excitations 
outside vortex cores. 

Collisions between the quasiparticles can cause some spreading of the incoming beam. 
However, the spreading can be made arbitrary small by lowering the density of the excita- 
tions, that is to say, by lowering the temperature. At low enough temperatures the mean 
free path exceeds the dimensions of the experimental cell and we can consider undamped 
excitations moving along straight paths until they hit a boundary or any potential barrier, 
particularly a barrier formed by a vortex. Andreev reflection of excitations thus gives the 
opportunity to probe flows in superfluid '^He at ultra-low temperatures. The most fruitful 
and promising application of Andreev scattering is thus superfluid turbulence in ^He-B in 
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the low temperature limit, that is to say for T/Tc < 0.4 K |17l . l29l . l43l . |44J . |45 |. 

Superfluid ^He-B is described by a macroscopic wave function, called the order parameter, 
with a well defined phase cj). The superfluid velocity is proportional to the gradient of 
the phase, 

V, = Av0, (2.1) 
Im 

where m is the mass of one ^He atom. Consequently, in contrast to classical fluids, superfluid 
motion is irrotational and vorticity exist only in the form of quantized vortices. Quantized 
vortices are line defects around which the phase (j) changes by 27r. The superfluid order 
parameter is distorted within the relatively narrow core of the vortex, and the superfluid 
flows around the core with speed which is inversally proportional to the distance from the 
vortex core. Since both the real and the imaginary parts of the order parameter are zero on 
the axis of a vortex, vortex lines can be considered as topological defects. Vortices cannot 
terminate in the middle of the flow, so they are either closed loops or extend to the walls. 

Superfluid turbulence consists of a tangle of quantized vortices. The complex flow fleld 
within the tangle acts as a potential barrier for quasiparticles, causing the Andreev reflection 
of a fraction of a beam of thermal excitations incident upon the tangle. The use of Andreev 
scattering as a visualization technique of ultra-low temperature turbulence requires to flnd 
out exactly what happens to a single quasiparticle which moves in the velocity fleld of a 
vortex, which is what we set out to do. 



III. EQUATIONS OF MOTION OF THERMAL EXCITATIONS 

Our flrst aim is to formulate, in the (x, ?/)-plane, the equations of motion of a single 
excitation moving in the velocity fleld of a single straight vortex which we assume to be 
fixed and aligned along the 2;-axis. We are thus concerned with a two-dimensional problem 
only. The quantities (here and below the numerical values of the quantities are taken at the 



bar pressure 46(]) which are necessary to describe the motion of the excitation are the Fermi 
velocity, vp ~ 5.48 x 10^ cm/s, the Fermi momentum, pp = Tifvp ~ 8.28 x lO"^*' g cm/s, 
and the Fermi energy, ep = pp'^ /{2m*) ^ 2.27 x 10~^^ erg. The quantity 



e. 



p2 



^ 2m* 



~ep (3.1) 
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is the "kinetic" energy of the excitation measured with respect to the Fermi energy, ep-, 
where m* 3.01 x m = 1.51 x 10~^^ g is the effective mass of the excitation, and p 
the momentum, p = |p|. Let Aq be the magnitude of the superfiuid energy gap. Near 
the vortex axis, at radial distances r smaller than the zero-temperature coherence length 
^0 = hvp/iT'^o ~ 0.8 X 10~^ cm, the energy gap falls to zero and can be approximated 



by A(r) ^ Aotanh(r/^o) |47|, |48|]. Since we are mainly concerned with what happens for 
r ^ C,Q, we neglect the spatial dependence of the energy gap and assume the constant value 
Ao = l.lQksTc ~ 2.43 x lO^^^ erg. 

Using polar coordinates (r, 0) in the (x, y) plane , the velocity field of a superfiuid vortex 
set along the 2;— axis is 

V. = ^e^, (3.2) 

h Tlh Q o , , . 

K = = — = 0.662 X 10-^ cmVs 3.3 

2m m 



where 



is the quantum of circulation, and is the unit vector in the azimuthal direction on the 
(x, ?/)-plane. 

In the presence of the vortex, the energy of the excitation becomes 



E = -^6^2 + A2 + p . V,. (3.4) 

In writing Eq. (13. 4p . the spatial variation of the order parameter is not taken into account 
for the sake of simplicity. We also assume that the interaction term p ■ varies on a spatial 
scale which is larger than ^o, and that the excitation can be considered a compact object 
of momentum p = p(t), position r = r(t), and energy E = E{p, r). This gives us the 



opportunity to use the method developed in Ref.[33], and consider Eq. (13.41) as an effective 
Hamiltonian, for which the equations of motion are 

dv dE{p, r) ep p 



dt dp A/e„2 + A^ rri* 



+ v„ (3.5) 



dp _ dE{p, r) _ ^ ) (3g) 

dt dr dr 

Eqs. (13. 5p and (13.61) have one immediate integral of motion, the energy: 

E{p, r) = E = constant. (3.7) 
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Eq. (13. 5p represents the group velocity of the excitation in the velocity field of the vortex. 
Excitations such that ep > are called quasiparticles, and excitations such that < are 
called quasiholes. The right-hand-side of Eq. (13.61) is thus the force acting on the excitation. 



IV. PROPAGATION OF EXCITATION IN THE VELOCITY FIELD OF A VOR- 
TEX 

We want to determine the trajectory of an excitation which moves in the two-dimensional 
velocity field of the vortex. It is convenient to rewrite the Hamiltonian, Eq. (13. 4p . and 
Hamilton's equations (13.51) and (13. 6p in polar cordinates (r, 0). We notice that the system 
consisting of a single excitation and a single vortex has a second integral of motion: the 
component of the angular momentum in the ^-direction, perpendicular to the plane of mo- 
tion, {x, y). Consequently, we can introduce two pairs of canonically conjugated variables, 
{pr'-i f) and (J = P(^r; 0), where Pr and p^j^ are the radial and azimuthal components of p, 
and J is the angular momentum. Since J is constant, it is convenient to write it in the form 
J = PfPo, thereby defining the constant po a particular trajectory. Under special initial 
conditions, as we shall see, po can be interpreted as the impact parameter. 



Eqs. (133!), dsn) and ([33]) become 




(4.1) 



(4.2) 



2m* 2m*r'^ 



dr _ ep Pr 



(4.3) 



dt ^ep2 + Al m* ' 




dt ^/cp^ + Aq 'm*r'^ 2nr'^ ' 



(4.4) 



By setting dE/dt = and using Eq. (14. 3 p we find 



^ dt 



dcp 



PfPo 



Pr K 



(4.5) 



m* irr 



3 ' 



and from Eq. (14.21) we have 




(4.6) 
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Eqs. f l4.ip - fl4.6p form a closed set which allows us to determine the trajectory of the 
excitation. 

It is apparent from Eq. fl4.3p that a quasiparticle incident upon the vortex has ep > 
and Pr < 0, whereas a quasiparticle moving away from the vortex has > and Pr > 0. 
Vice- versa, a quasihole incident upon the vortex is characterized by < and Pr > 0, 
whereas a quasihole moving away from the vortex has ep < and pr < 0. 

Later we shall consider a quasiparticle which leaves a point of the wall of the cylindrical 
experimental cell; this quasiparticle is initially characterized hj r = R (where R is the radius 
of the cell), p = Pf and pr < 0. The axis of the vortex will still be at the centre of the 
coordinate system. In such case the quasiparticle with initial momentum directed along the 
X-axis will feel the effective pairing potential Aejj ^ Aq — pFyn/{2TTr) (Fig. [T]). 

It is obvious from Eq. (14.30 that unless po is exactly zero (J = 0), the radial velocity 
of the excitation will eventually vanish. This may happen either because Pr = (classical 
turning point) or because ep = (Andreev turning point). 

It can be seen from Eqs. (14. ip and (14. 6 p that the classical turning point is reached first 
when 

E>Ao+Pf^^aJi + ^-^) (4.7) 

(here and in the equations below the numerical factor 3 is introduced by the ratio between 
the effective mass of quasiparticle and the bare mass of a ^He atom: m*/m ^ 3) in which 
case a quasiparticle with this energy follows a trajectory which is of the "normal" type: the 
quasiparticle retains its "particle" nature and moves past the vortex, across the experimental 
cell to the wall on the opposite side. On the contrary, a quasiparticle with energy E such 
that 

Ao<E<Ao + PF-^^Ao(l + ^) (4.8) 

reaches the Andreev turning point first, undergoes Andreev reflection, and returns to a point 
near its starting point after changing its nature and becoming a quasihole. 

Of these two cases, our concern is the case of Andreev reflection. We first determine the 
locus of Andreev turning points, defined by the minimum radial distance from the vortex 
core: 

'-"^ - [2n {E - A,) J -{ 2 (E-Ao)) ' ^^''^ 
Consider a quasiparticle which has reached r = rmin- At this point the radial velocity 
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r vanishes, but the excitation does not stop. It has still a nonzero azimuthal velocity, rcf). 
Thereafter the excitation propagates as a quasihole (characterized by a negative value of e^) . 

In order to calculate the trajectory of a reflected quasiparticle it is convenient to simplify 
the governing equations of motion using the fact that at the ultra-low temperatures which 
interest us, T <^ Tc, most quasiparticles have energies ep ^ Aq. We can then make the 
following approximation: 



This spectrum is similar to Landau's spectrum of excitations in superfiuid He II near the 
roton minimum {p = po), E ^ Aq + {p — PoY/ (27/1^) (where is the effective roton mass). 



which was used to calculate the mutual friction force 49|]; note that in Eq. f l4.10p the role 
of the roton mass is played by the ratio Ao/vp. 

Using Eqs. fl4.10p . (14 .Sp and (I4.4p . and the smallness of the ratios ep/ep and Ao/ep, we 
obtain 

_ m*Ao _ m* Tmin r^dr {A^^\ 

~ eppr "~ PF {3ir^oPoy/Hr^ - rlJ^r' - Ply/' ' ^ 

^0 ± ^^''"''^ ^ dr U 12) 



r\\ - pg/r2)i/2 2(3vreoPo)i/^(r2 - r^,jV2(r2 - p2)i/2^ 
where h = h/pp, the sign plus is used for quasiparticles and the sign minus for quasiholes. 

From Eqs. (14.111) and (14.121) we obtain the Andreev return time r of the excitation (the 
time it takes to travel from the radial distance R to the Andreev reflection point and back) 
and the Andreev reflection angle A0: 

~ fn-^ „ M /9 / 7Z3. I2 VUTTj?. „2M 19. ' (4-13) 



VF (3<oPo)^/^ (r^ - rLJ^/^(r^ - pg) V 

The evaluation of these elliptic integrals is shown in the Appendix. We obtain 

, _ .r^^n 1 ( {R^- rLnmn' - Plf" , , .4l5^ 

VP {^^i,P,yl'\ R ' ^ ^ 

which becomes, assuming R^ r and R^ po, 

Rr^in R f 2Ao \ R 2 Aq 



t^F(3<oPo)i/2 vf\E-AoJ VF ep ' ^' ^ 
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We conclude that the Andreev return time is longer if the excitation's energy is lower. 
Similarly, assuming po/R <ti 1 and rmin/R 1, the Andreev reflection angle is 

To apply these results we assume that the initial momentum of the quasiparticle is di- 
rected along the x-axis, and that the angular momentum J = ~pyo = PfPo- From Eq. (13.1 P 
it follows that the momentum p = ppi^ + 2m*ep/ppY^'^ and, in the ultra-low temperature 
limit, {p — Pf)/pf < 10~^. For we have = po(l + 2m*ep/pF)~^^'^ ~ Po- In this 
case Po becomes the impact parameter (Fig. [2]), and Eq. (I4.17P shows that quasiparticles 
with smaller impact parameter (hence smaller angular momentum) are Andreev reflected 
by smaller angles. 

As it is seen from Eq. (14. 9p . the Andreev radius depends strongly on the initial energy of 
the excitation: 

r-mm = (3<oPo)'/'— • (4.18) 

The same arguments apply to the critical value poc defined as a maximum value of po 
which causes the Andreev reflection of quasiparticles with the given initial energy ep. To 
calculate poc, we assume that at the starting point of the trajectory the quasiparticle has 
coordinates {R, 0o); where 0o = arcsin(?/o/-R) ~ — arcsin(poc/-R); the coordinates of the 
Andreev reflection point in this case should be (r^m, ^7r/2). Thus the difference between 
the reflection angle and the starting angle is 

A(f) = -- + arcsin (^) . (4.19) 



2 \ R 

This difference can also be calculated from Eq. (14. 12p where the second term (of the order 
of h/pp) in the integrand can be neglected. We obtain 

A0 = — arcsin ( | -)_ arcsin ) . (4.20) 

By comparing Eqs. fl4.19p and (I4.20p we find 

Poc ~ 3<o (^^^ . (4.21) 

In the typical low temperatures experiments ksT / /S.q ^ 0.1, and, for quasiparticles with 
initial energy ep ~ ksT, we find r^j^ ^ 10(37r^oPo)^/^ and Poc ~ lO^^o, while the same 
quantities for the quasiparticles with Cp ~ (AoA;bT)^/^ are r^m ~ 3(37r^oPo)^''^ and poc ~ 



V. HEAT TRANSPORT THROUGH THE VELOCITY FIELD OF A VORTEX 



In the experimental studies of superfluid turbulence in ^He-B at the ultra-low temper- 
atures the vortex tangle is studied by detecting the fraction of quasiparticles which are 
Andreev reflected by the vortices and measuring the heat which is transported by the quasi- 
particles. Using the results of previous Sections, it is straightforward to calculate the fraction 
of energy (or heat) transmitted across the velocity field of a vortex. Once we know this frac- 
tion, we shall generalize it to a system of many vortices. 

In Section IV it was explained that the quasiparticles characterized by the particular 
impact parameter po a-re Andreev reflected by a vortex if their energies satisfy the condition 
Ao < -B < Ao(l-|- 37r^o/2po)- If this condition is not satisfied, the quasiparticles pass freely 
across the vortex velocity field. If the system which consists of the vortex and quasiparticles 
is in thermal equilibrium, there is no preferred direction around the vortex. Incident and 
transmitted fluxes at one side of the vortex are canceled by the fluxes in the opposite 
direction, and no net flow of energy exists when the temperature everywhere around the 
vortex has the same value. 

A net flux of quasiparticles and of energy results only if there is some small temperature 
difference, 5T <^T between the two sides. If this is the case, the heat carried by the incident 
quasiparticles is described by the expression: 

SQinc = N{E)vg{E)E^^ STdE, (5.1) 

where 



Here 



= ^ (5-3) 



is the density of states at the Fermi energy with corresponding Fermi momentum pp. The 
group velocity of a Bogolubov quasiparticle Vg is given by the expression: 



^. = 1^^- ^ VF, (5.4) 



and f{E) is the Fermi distribution function, which, at ultra-low temperatures, is transformed 
into the Boltzman distribution, and describes the mean occupation number of a state with 
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energy E: 

f{E) = e~^. (5.5) 

The thermal flux of quasiparticles incident on the vortex velocity field per unit length 
per unit time is obtained with the help of Eqs. (15. 2p . (15. 4p and (15.51) : one has 

roc ^ 

^Q^-- = ^P''Pk^ ^V^rfE^iV^t;^A2— e-^. (5.6) 

If there is a plane current of quasiparticles with transverse cross section Rq, then the total 
heat current incident on the vortex per unit time will be: 

5T aq 

Qinc = 2Ro5Qinc = '2RoNFVF^lY^ (5-7) 

We assume that, in the (x, |/)-plane orthogonal to the straight vortex line, the polar- 
ity of the vortex located at (0, 0) is positive and consider quasiparticles incoming in the 
positive X— direction. As discussed earlier, in this case the upper half-plane will be abso- 
lutely transparent for quasiparticles so that the heat transferred by quasiparticles through 
this half-plane will meet no resistance. The lower half-plane of vortex flow field will reflect 
a fraction of quasiparticles and induce some thermal resistance. A quasiparticle with the 
impact parameter po is transmitted through the vortex velocity field if it carries the energy 
E > Ao(l + |7r.^o/Po)) in which case the heat which is transferred per unit time by such a 
quasiparticle can be calculated as 

/■°° df(E) \ ( 37r£n\ ^0 

<5Q(p)=/ Ar(E)t;,(£;)E-LLi5Trfi?^g_-- i + e-^^. (5.8) 

Notice that estimating the ratio .^o/Po we kept only the linear term. 

The total amount of energy transferred through the vortex by quasiparticles originated 
within the interval —Rq ^ y < Rq is: 



Qtr 



2 

The integral in Eq. (l5.9p can be estimated as 



Ro Jo \ Po J 



(5.9) 



^ Roe ''bt . (5.10) 
Thus the fraction of heat which is transferred through the velocity field of the vortex is 

1 / _^3!L€o\ 

Sft^ = - il + e "B^ 2flo . (5.11) 



2 
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In experiments at ultra-low temperatures we have Ao/kBT ~ 10, so that the cross-section 
of the thermal flux is Ro ~ lO^o and approximately 52% of the total heat is transferred 
through the vortex. If the heat current has the cross-section ~ lO^^o^ the fraction of the 
transferred heat is approximately 0.82. Therefore the reflection of the heat flux takes place 
only in the close vicinity of the vortex core. 



VI. ANDREEV REFLECTION IN A VORTEX GAS. 

To apply our result to experiments, we consider for simplicity a system of random parallel- 
antiparallel vortices (i.e. a system of vortex points in the (x, ?/)-plane; such a system is known 
as the Onsager vortex gas). This vortex system is penetrated by a quasiparticle current 
created by a temperature difference 6T. It is convenient to introduce the effective radius 
Ro of each vortex as the half of the mean intervortex distance, i.e. Rq = i/2. We divide 
the vortex configuration in parallel layers of width i each perpendicular to the quasiparticle 
current. Clearly, the transmittability of each layer is equal to the transmittability of a vortex 
within a region of radius i/2. Thus the fraction of heat transmitted by each layer is 

5f,, = -{l + e . (6.1) 



2 

If we assume now that vortices are well separated and that their velocity fields do not 
overlap significantly, we obtain the conditions 



Eq. (16.11) becomes 



2\ ksT i J kBT 2i ^ ' 

Driven by the temperature difference, the heat flux Qq reduces, after penetrating the first 
layer, to Qi = QoSftr', after penetrating the second layer, it becomes Q2 = QiSftr- Hence, 
after penetrating the last nth layer, we obtain Qn = Qn-i ^ftr- Thus we have 

Qn = Qn-lSftr = ... = QoSfTr- (6-4) 

We conclude that the fraction of heat which is transferred through the system of vortices is 

ftr = {5Ur. (6.5) 
12 



If the total vorticity is confined within a region of size 5", the number of layers, n can be 
estimated ~ S/i. From Eq. (16.51) we obtain: 



fi 



tr 



Finally we obtain the intervortex distance: 



Ao Svr^o'S' 



(6.6) 



(6.7) 



ksT 2 In ftr ^ 

The quantities S (the size of the vortex system) and ftr (the fraction of reflected quasi- 
particles) in E q. (6.7ll can be observed experimentally. From the available description of 
291], the maximum transmitted fraction of quasiparticle current is ftr ~ 0.75 



one experiment 



and the spatial extension of the vorticity is S* ~ 2 ■ 10^^ cm. Since the zero temperature 
coherence length is ~ 0.8 x 10"^ cm, we conclude that in the case where Aq/Z^bT ~ 10 the 
average intervortex distance is £ ~ 1.62 ■ 10^^ cm, which is in good agreement with existing 



estimates 



43|. 



VII. CONCLUSIONS 



In conclusion, starting from Hamilton's equations, we have calculated the trajectories of 
quasiparticles which move in the velocity field of a quantized vortex in ^He-B and determined 
the Andreev reflection point. Generalizing the result to a disordered system of many vortices, 
we have determined the precise location of turning point and showed how to recover the 
typical intervortex spacing in the turbulent '^He-B. Our result is in good agreement with less 
rigorous estimates. 

Future work will investigate Andreev reflection of quasiparticles by a system of moving 
vortices. We shall also study how the Andreev reflection technique can be used to visualize 
vortex structures (e.g. coherent bundles of vortices) and determine turbulent fluctuations 
and turbulence statistics. 
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APPENDIX A 



The Andreev return time, r, and the Andreev reflection angle, A0, are defined by for- 
mulae (14. 131) and fl4.14p . where < po < ^min < r < R. To evaluate these formulae we use 
the following integrals: 

dr ^JB? — 



h 



R 



(A.l) 



where 

G = K 
and 



\ ' mm / \ 



arcsm 



f^min Po 



B ' r 



-E 



' mm / \ 



1 J?0 

2' r„ 



arcsm 



R r. 



min 



(A.2) 



dr 



where K, F and E are elliptic integrals, defined as 

F{k 0) = 



K (-^) - F ( 

\ ' mtn / \ 



arcsm 



^min Po 



B r 

^ ' m.ir. 



(A.3) 



V 1 - A;2 sin^ (0) 



K{k) = F A:) , 
E{k, = ^l-F sin2(^), 



(A.4) 

(A.5) 
(A.6) 



with 6 = arcsin {rmml R) and k = pq/R. 

For k"^ < 1 the elliptic integrals flA.4p . flA.Sp and flA.6p are represented by the series 



20 k^ 
F{k, 9) = —K{k) - sine cos + 

TT 4 



K{k) = l + \k^ + ..., 

26 k"^ 
E{k, 6) = —E{k) + sin 6 cos d— + 

E(k^^^^E,k)^l-lk^,... 



(A.7) 

(A.8) 

(A.9) 
(A.IO) 



using which we obtain 



2r 



T 



2 

min 



4 Vr, 



(A.ll) 
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Assuming R 3> Tmin, Po and po < ^min we have 

2/? r 

^ ^ z-n, ; (A.12) 

Similarly, 
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FIG. 1: (Color online) The dimensionless effective potential D = Ae//(r/^o)/^o seen by quasi- 
particles with momentum parallel to the x-axis and moving from x = — oo. The dimensionless 
coordinates x and y are in units of ^o- 



18 




FIG. 2: (Color online) Schematic trajectory of the quasiparticle which starts at position A, is 
Andreev-reflected by the vortex (at the origin) at position B (where it becomes a quasihole), then 
traces its way back with a small Andreev angle A(f) (not to scale). 
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